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ABSTRACT 

The  effects  that  regional  size  and  various  predictive 
variables  have  on  predicting  nonindustrial  private  timber 
harvest  in  Montana  are  presented.  Multivariable  models 
predicted  harvest  substantially  better  than  single 
variable,  price-based  models.  Multicounty  area  models 
were  not  consistently  better  predictors  of  harvest  than 
the  State-level  model.  Considerable  losses  in  predictive 
ability  relative  to  area-specific  models  were  found  when, 
for  consistency  purposes,  a  common  model  was  used  to 
predict  harvest  for  each  area. 


KEYWORDS:  timber  harvest,  nonindustrial  private 
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When  one  thinks  of  timber  harvest  on  nonindustrial 
private  forest  land  (NIPF),  visions  of  the  eastern  part  of 
the  country  come  to  mind— the  Northeast,  Midwest,  and 
South— certainly  not  the  Intermountain  West  or 
Montana.  These  lands  are  dominated  by  public  forests 
and  NIPF  harvest  is  of  little  concern.  Right?  Well,  only 
partially  so.  While  a  decade  and  a  half  ago  the  NIPF 
sector  produced  about  2  percent  of  Montana's  timber 
harvest,  it  now  produces  one-fifth  of  the  State's  harvest 
or  about  200  million  board  feet  annually  (unpublished 
data,  Northern  Region,  Forest  Service,  and  Montana 
Division  of  Forestry).  Montana's  increased  NIPF  harvest 
reflects  a  more  general  realinement  of  timber  supply 
sources  in  the  Intermountain  region— emphasize  private, 
deemphasize  public. 

With  the  increasingly  important  role  played  by  the 
NIPF  base  has  come  increasing  interest  in  better 
understanding  and  predicting  timber  harvests  from  these 
lands.  Concerns  historically  shown  by  State  forestry 
organizations  are  now  being  expressed  by  the  Forest 
Service,  U.S.  Department  of  Agriculture,  where  the  need 
for  forest-  and  regional-level  planning  compels  more  ex- 
plicit attention  to  NIPF  timber  production. 

Unlike  the  long  history  of  NIPF-oriented  investiga- 
tions found  elsewhere,  the  Intermountain  West  and  Mon- 
tana have  few  or  no  such  traditions.  Recently  available 
timber  harvest  and  supply  analyses  pertaining  to  these 
areas  are  not  well-suited  to  needed  evaluations  (see 


Research  forester  and  economics  assistant,  respectively,  at  the 
Intermountain  Station's  Forestry  Sciences  Laboratory,  Missoula,  Mont. 


USDA  Forest  Service  1982  and  Adams  and  Haynes 
1980).  Because  these  analyses  have  been  directed  toward 
a  national  timber  assessment  with  eight  major  regions, 
relationships  for  smaller  geographical  areas  were  not 
developed.  Additionally,  the  supply  data  for  the  Rocky 
Mountain  region  applied  to  the  total  private  sector,  not 
distinguishing  forest  industry  from  NIPF. 

The  research  reported  here  was  conducted  in  response 
to  a  continuing  need  to  predict  timber  harvest  for  all 
producing  sectors  of  the  State  and  sub-State  levels.  The 
research  hypothesis  asserted  that  substantial  im- 
provements could  be  made  by  estimating  NIPF  timber 
harvest  separately  from  forest  industry,  by  adopting 
smaller  geographical  regions,  and  by  expanding  the 
range  of  predictor  variables  available.  While  appearing 
similar  to  models  of  economic  supply,  models  developed 
were  formulated  differently  to  predict  timber  harvest 
only. 

METHODS 

Montana  was  chosen  for  study  because  of  its  impor- 
tance in  Intermountain  States  timber  production  and 
readily  available  data.  Figure  1  shows  the  boundaries  of 
the  four  geographical  regions  used  in  this  study.  Based 
on  available  timber  production  and  marketing 
characteristics,  two  small  multicounty  areas  were 
identified— northwestern  and  southwestern  Montana. 
Either  contains  a  land  base  greater  than  Maryland  or 
Vermont.  These  two  areas  also  were  aggregated  into  a 
western  Montana  region,  with  a  land  area  about  the  size 
of  West  Virginia.  The  final  region  consisted  of  the  total 
State,  the  fourth  largest  in  the  Nation. 

Selection  of  the  potential  independent  variables  such 
as  price  was  based  on  economic  theory,  while  others 
such  as  ownership  size  were  based  on  previous  NIPF 
studies  that  identified  factors  affecting  landowner  timber 
harvest  behavior.  Five  broad  categories  were  identified: 
(a)  characteristics  of  the  NIPF  land  base,  (b)  economic 
conditions  prevailing  in  the  area,  (c)  price  received  for 
the  stumpage,  (d)  timber  harvest  behavior  of  other  stum- 
page  producers,  and  (e)  a  miscellaneous  category. 

While  theory  and  previous  studies  suggested  these 
categories,  neither  adequately  specified  the  actual 
variables  to  be  used.  Consequently,  we  selected  several 
candidate  variables  for  each  category  that  had  the 
potential  of  being  useful  in  prediction.  Further,  we  had  a 
choice  of  expressing  many  of  these  candidate  variables 
in  several  ways;  for  instance,  dollar-measured  candidate 
variables  could  be  expressed  in  nominal  or  constant 
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Figure  1.— Montana  and  three  multicounty  regions  used  to  assess  prediction  models  in  this  study. 


(real)  dollars  and  most  variables  could  be  expressed  in 
current  year  or  previous  year  (lagged  1  year)  levels.  Con- 
stant dollar  conversions  used  the  Gross  National  Prod- 
uct Implicit  Price  Deflator,  1972  base.  We  ultimately 
treated  variables  in  one  of  three  ways: 
Type  Treatments 

I  Nominal  dollars;  current  year 

II  Nominal  and  real  dollars; 

current  and  previous  years 
III  Nominal  dollars;  current  and 

previous  years 
Table  1  summarizes  the  categories  of  variables,  types  of 
candidate  variables,  and  a  treatment  code  (I,  II,  or  III) 
for  each.  In  fact,  we  used  even  more  candidate  variables 
than  are  indicated.  For  example,  the  per  capita  income 
variable  shown  was  really  represented  by  three  variant 
forms:  per  capita  personal  income,  per  capita  disposable 
personal  income,  and  per  capita  nonfarm  income,  each 
treated  in  a  Type  II  manner.  In  total,  we  quantified 
several  dozen  analysis  variables,  combinations  of  can- 
didate variables,  variant  forms,  and  methods  of 
treatment. 


Data  series  were  developed  for  the  dependent  variable, 
NIPF  timber  harvest  (MMBF,  Scribner) ,  and  for  each  of 
the  potential  independent  variables,  for  each  designated 
region,  for  the  16-year  period  1965  through  1980. 
County-level  or  quasicounty-level  data  were  generally  ob- 
tained and  aggregated  into  study  regions.  Savings  rates, 
interest  rates,  and  two  price  indexes  are  national  in 
scope.  A  wide  variety  of  data  sources  were  used. 

Data  were  fit  to  multiple  linear  regression  models  us- 
ing ordinary  least  squares  and  the  BMDP  stepwise 
regression  computer  procedure  (Dixon  1981).  In  linear 
regression,  available  data  are  used  to  estimate  coeffi- 
cients for  the  model's  independent  variables  so  as  to  ac- 
count for  as  much  variability  in  the  dependent  variable 
as  possible.  Region- specific  data  for  all  available  years 
were  used  to  estimate  each  model.  Development  of  the 
"best"  model  proceeded  until  either  no  new  variables 
were  statistically  significant  (a  —  0.1)  or  when  the  three 
most  significant  variables  had  been  identified  (see 
Draper  and  Smith  1981). 

For  the  models  reported  in  this  paper,  the  un- 
restricted, automatic  stepping  process  of  the  computer 


Table  1  .—Categories  of  independent  variables,  candidate  variables,  and  treatment  code  (I,  H,  or  III)  used  to  predict  NIPF 
timber  harvest 


NIPF  characteristics 

Economic  conditions 

Price 
received 

Other  harvests 

Miscellaneous 

Acres  per  owner,  I 

Per  capita  income,  II 

$/MBF,  II 

National  Forest,  III 

Lumber  output,  III 

Volume  per  acre,  I 

Interest  rate,  national,  III 

Forest  industry,  III 

Residual,  III 

Growth  per  acre,  I 
Timber  size,  I 

Savings  rate,  national,  III 
Unemployment  rate,  III 

Uncut  volume  under 
contract,  III 

procedure  was  frequently  interrupted  to  preclude  selec- 
tion of  "redundant"  variables.  Many  candidate  variables 
were  handled  with  variant  forms  and  methods  of  treat- 
ment because  of  uncertainty  as  to  which  form  or  treat- 
ment would  prove  most  useful  statistically.  When  one  of 
these  variables  was  selected  by  the  computer  program 
for  inclusion  into  a  model,  the  related  variables  were 
then  irrelevant,  superfluous  to  the  analysis,  and  exclud- 
ed from  subsequent  consideration.  For  example,  if  the 
per  capita  personal  income  variable  (real  dollars,  current 
year)  was  selected  for  inclusion  in  a  model,  not  only  were 
all  other  methods  of  treatment  for  that  variable  (for 
example,  nominal  dollars,  previous  year)  excluded  from 
future  consideration,  but  the  alternative  variant  forms 
(such  as,  per  capita  disposable  personal  income  and  per 
capita  nonfarm  income)  were  excluded  as  well.  Once  the 
"best"  models  were  identified,  a  series  of  "all  possible 
combinations"  of  variables  regression  runs  was  made  to 
determine  if  a  preferable  model  existed  but  was  missed 
by  the  stepwise  procedure. 

All  "best"  models  were  tested  for  statistical 
significance  and  consistency  with  underlying  assump- 
tions. Statistical  significance  of  regression  coefficients 
was  determined  with  a  two-tailed  test.  That  procedure, 
which  ignores  the  sign  of  the  coefficient,  was  judged  ap- 
propriate for  prediction-motivated  modeling.  Subsequent 
tests  for  multicolinearity,  heteroscedasticity,  and 
autocorrelation  either  revealed  no  problem,  or  problems 
were  resolved;  one  autocorrelation  test  was  inconclusive 
(see  Koutsoyiannis  1977).  Results  presented  use  the 
unadjusted  coefficient  of  multiple  determination  (R2) 
rather  than  the  adjusted  coefficient  because  of  its 
understandability  while  faithfully  reflecting  the  adjusted 
coefficient  in  these  analyses. 

RESULTS 

Predictive  equations  resulting  from  our  modeling  ef- 
forts seemed  quite  satisfactory,  accounting  for  three- 
fourths  to  almost  90  percent  of  the  variation  in  NIPF 
timber  harvest.  Equations  for  these,  "best"  models  are 
shown  below  for  Montana  (MT),  western  Montana  (WM), 
southwestern  Montana  (SW),  and  northwestern  Montana 
(NW)  regions  together  with  t-statistics,  level  of 
significance  (a  =  0.01,  0.05,  or  0.10),  and  the  percentage 
of  variation  in  harvest  explained  by  the  model,  the  R2. 

H^T  =  258.89  +  4.75(P)t  -  0.16(Lmbr)t_j  .  .  .  R2  =  80.7 

(6.72,  .01)     (-2.95,  .05) 

H™=   -27.99  +  2.70(P)t  +  8.47(8),..!  .  .  .  R2  =  83.7 

(8.15,  .01)  (-2.60,  .05) 

nt    =     95.71  +  0.76(P)t  -  1.27(Res)t  +  2.4(Unemp)t.  .  .  R2  =  89.3 
(5.22,  .01)  (-5.29,  .01)       (1.86,  .10) 

uNW 

Ht    =     -8.37  +  1.30(P)t  +  0.80(S)2-i  .  .  .  R2  =  76.5 

(5.59,  .01)    (4.48,  .01) 

where: 

H  =  NIPF  timber  harvest;  region  (MMBF,  Scribner). 
P        =  National  Forest  stumpage  price,  cut;  region; 

$/MBF  (1972  $). 
Lmbr  =  lumber  output;  region  (MMBF). 
S         =  savings  rate;  nation  (percent  =  savings  as 

percent  of  personal  disposable  income). 


Res      =  residual;  region  (percent  =  National  Forest 
harvest  +  industry  harvest  as  percent  of 
lumber  output). 

Unemp=  unemployment  rate;  region  (percent  =  un- 
employed workers  as  percent  of  labor  force). 

t         =  year. 

Given  the  large  number  of  candidate  independent 
variables,  the  recurrence  of  certain  variables  in  the 
"best"  models  is  comforting,  if  not  somewhat  surprising. 
The  real  price  received  for  National  Forest  stumpage  cut 
was  clearly  the  single  best  predictor  variable  in  all  four 
models.  We  do  not  contend  that  National  Forest  stum- 
page prices  cause  NIPF  owners  to  cut  timber;  clearly 
they  do  not.  Rather,  because  no  data  are  available  on  ac- 
tual prices  paid  to  NIPF  owners,  any  available  price 
data  can  only  represent  an  index  or  proxy  for  the  actual 
stumpage  price  received  by  the  landowners.  National 
Forest  stumpage  cut  prices  and  NIPF  stumpage  prices 
probably  move  together. 

In  addition  to  price,  the  national  savings  rate  (S)  was 
an  important  predictor  variable  in  the  western  and 
northwestern  Montana  models.  In  both  cases,  the 
"previous  year"  form  of  the  savings  rate  variable  was 
used  and  the  sign  of  the  coefficient  was  positive.  The 
residual  variable  (Res)  was  also  an  important  predictor 
in  southwestern  and  Montana  models.  In  the  "current 
year"  form,  it  had  a  negative  coefficient.  Loosely  inter- 
preted, this  variable  represents  (on  a  percentage  basis) 
the  proportion  of  the  region's  lumber  output  accounted 
for  by  timber  harvested  from  industry  and  National 
Forest  lands.  The  higher  the  proportion,  the  lower  the 
NIPF  harvest.  It  portrays  NIPF  harvest  as  filling  an 
unmet  or  residual  need  for  timber.  The  "previous  year" 
form  of  the  lumber  output  variable  (Lmbr)  was  signifi- 
cant in  the  Montana  model.  Finally,  unemployment  rate 
was  important  in  the  southwestern  Montana  model,  be- 
ing positively  related  to  NIPF  harvest.  The  higher  the 
rate  of  unemployment,  the  higher  the  harvest. 

All  "best"  models  accounted  for  a  substantial  portion 
of  the  variation  in  NIPF  timber  harvest,  as  reflected  by 
the  relatively  large  R2  values.  The  larger  this  value,  the 
more  variation  explained  by  the  model,  and  the  more 
useful  the  model  is  in  prediction.  Table  2  shows  a  com- 
parison of  the  actual  NIPF  harvest  and  the  harvest 
predicted  by  the  "best"  model,  for  each  region.  The 
"best"  models  vary  in  prediction  ability.  The  north- 
western, southwestern,  and  western  Montana  models 
had  the  smallest  average,  annual  percentage  difference 
between  observed  and  predicted,  each  averaging  just 
under  10  percent  over  the  16-year  period.  The  Montana 
model  had  the  largest  differences,  averaging  about  13 
percent  per  year  over  the  period.  These  larger  differences 
are,  no  doubt,  related  to  the  volatility  of  eastern  Mon- 
tana harvest. 

Derivation  of  predicted  harvest  levels  can  be  il- 
lustrated by  the  southwest  model  for  1980.  Given  the 
1980  level  for  real  price  (P  =  $20.87/MBF),  residual  (Res 
=  77.0  percent),  and  unemployment  (Unemp  =  7.5  per- 
cent), the  level  of  harvest  is  predicted  as  follows: 

HfWQ    =  95.71  +  0.76(20.87)  -  1.27(77.0)  +  2.4(7.5) 
=  31.4  MMBF 
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Table  2.—  Comparison  of  actual  NIPF  timber  harvest  and  that  predicted  by  "best"  models,  by  region, 
1965-80 


Northwestern  Southwestern  Western  Montana 

Year      Predicted      Actual      Predicted      Actual      Predicted      Actual      Predicted  Actual 


1965 

43.6 

59.0 

1966 

41.9 

36.7 

1967 

44.7 

42.8 

1968 

61.2 

49.3 

1969 

70.9 

63.9 

1970 

49.4 

49.0 

1971 

61.7 

72.5 

1972 

82.1 

92.1 

1973 

69.6 

75.2 

1974 

86.4 

86.7 

1975 

64.6 

53.3 

1976 

70.5 

67.2 

1977 

61.3 

61.3 

1978 

63.2 

65.4 

1979 

65.9 

62.9 

1980 

42.6 

42.0 

Million  board  feet- 


19.9 

18.7 

64.3 

13.3 

12.1 

63.1 

18.5 

19.1 

65.1 

23.8 

21.8 

84.1 

35.9 

45.6 

117.1 

43.3 

39.3 

84.9 

24.6 

27.1 

86.5 

34.3 

30.7 

122.1 

56.8 

53.6 

118.8 

40.1 

46.5 

134.9 

37.8 

30.7 

95.2 

47.9 

49.5 

105.8 

56.4 

64.6 

110.1 

49.9 

51.5 

114.5 

57.7 

44.9 

129.9 

31.4 

31.2 

74.6 

77.7 

112.3 

124.4 

48.8 

108.5 

105.1 

61.8 

99.8 

98.4 

71.1 

121.5 

99.9 

109.5 

160.3 

151.2 

88.3 

138.9 

140.0 

99.6 

140.3 

180.1 

122.8 

173.1 

134.9 

128.9 

200.6 

215.7 

133.2 

184.9 

229.7 

84.0 

176.2 

148.7 

116.7 

214.7 

231.1 

125.9 

222.0 

218.2 

116.9 

227.6 

222.2 

112.9 

238.3 

227.9 

73.2 

157.4 

148.6 

Estimates  for  other  regions  and  years  are  similarly 
derived.  The  estimate  for  the  1981  harvest,  a  year  out- 
side the  modeling  data  base,  in  the  southeastern  region 
is  calculated  as  follows: 

H?98i    =  95-71  +  0.76(21.85)  -  1.27(86.1)  +  2.4(9.2) 
=  25.0  MMBF 
The  actual  level  of  the  1981  harvest  for  that  region  was 
25.1  MMBF.  Differences  between  the  predicted  and  ac- 
tual harvests  for  the  northwestern  and  western  Montana 
regions  were  typical  of  the  16-year  period;  that  for  Mon- 
tana exceeded  the  typical  difference. 

One  purpose  of  our  work  was  to  see  if  improvements 
in  predicting  NIPF  harvest  could  be  made  by  expanding 
the  range  of  independent  variables  considered.  We  found 
that  although  addition  of  the  second  and  third  variables 


could  markedly  improve  prediction,  choice  of  the  first 
variable  was  critical.  Table  3  shows  a  comparison  of 
model  explanation  capability  (R2)  between  single  variable 
models  using  several  real  price-based  variables  and  the 
multiple-variable  "best"  models.  As  shown,  im- 
provements can  be  made  over  single,  price-based  variable 
models  by  including  more  explanatory  variables.  The 
largest  improvement,  a  90  percent  increase,  was  obtain- 
ed in  the  two-variable  northwestern  model,  relative  to 
the  best  price-based  model.  This  point  is  further  em- 
phasized when  the  multivariable  models  are  compared  to 
a  widely  used  price  variable,  National  Forest  stumpage 
(sold)  price. 

Substantial  prediction  improvements  can  be  made  by 
using  the  "correct"  price  variable.  Inspection  of  the 


Table  3. — Comparison  of  explanatory  ability  (R2)  for  alternative  "real"  price-based  models 
relative  to  "best"  multivariable  models,  by  region 


Alternative  price  measures/indexes1 


"Best"  model;2  variables  = 


Region 

WWPA3 

WPI4 

Mont5 

NFS6 

NFC7 

1 

2 

3 

Montana 

18.5 

52.0 

57.2 

30.7 

67.7 

67.7 

80.7 

Western 

0.3 

22.7 

49.8 

55.2 

75.2 

75.2 

83.7 

Southwestern 

7.0 

44.2 

49.8 

50.3 

55.5 

55.5 

86.2 

89.3 

Northwestern 

2.7 

2.5 

25.4 

16.8 

40.2 

40.2 

76.5 

1Equivalent  models  using  "current"  prices  were  substantially  inferior  to  "real"  price  models. 

independent  variables  significant  at  a  =  0.10  level,  or  better. 

3Western  Wood  Products  Association,  price  index,  dry  Douglas-fir  and  larch,  real. 

4Wholesale  Price  Index,  lumber  and  wood  products,  1972  base. 

5Montana  Division  of  Forestry  stumpage  price,  sold,  S/MBF,  real. 

6National  Forest  stumpage  price,  sold,  S/MBF,  real. 

7National  Forest  stumpage  price,  cut,  S/MBF,  real. 
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predictive  power  of  alternative  price  measures/indexes 
shows  wide  variation.  The  Western  Wood  Products 
Association  price  index  is  generally  the  worst,  and  the 
National  Forest  stumpage  (cut)  price  is  always  the  best. 
In  all  cases,  real  prices  were  found  to  be  far  more  useful 
in  prediction  than  their  nominal  price  counterparts.  All 
this  suggests  a  need  for  analysts  to  rather  carefully 
choose  type  of  price  variable  to  use. 

Another  purpose  served  in  this  study  was  to  deter- 
mine the  extent  to  which  improvements  in  predicting 
NIPF  harvest  could  be  secured  by  focusing  on  smaller 
geographical  areas.  Frankly,  we  expected  substantial  im- 
provements. But  this  was  not  the  case.  Table  3  shows 
mixed  results.  A  comparison  between  the  "best"  one-, 
two-,  or  three-  variable  models  associated  with  increas- 
ingly larger  regions  (NW  and  SW,  to  WM,  to  Montana) 
shows  no  consistent  pattern.  Models  for  the  smaller 
regions  were  better  or  worse  than  those  for  large  ones, 
depending  on  the  specific  regions  being  compared  and 
the  number  of  variables  in  the  models.  These  results  are 
suggestive  only,  definitive  conclusions  being  possible 
only  from  more  extensive  research  involving  replications. 
Nevertheless,  we  speculate  that  major  improvements  in 
predictability  were  secured  by  initially  focusing  on  a 
single  State  rather  than  a  large,  multi-State  region. 

Timber  harvest  for  western  Montana  can  be  predicted 
in  one  of  two  ways.  It  can  be  predicted  directly  by  use 
of  its  "best"  model  as  described  earlier.  It  can  also  be 
predicted  by  summing  the  predicted  harvests  for  the 
northwestern  and  southwestern  regions.  These  summa- 
tions provided  a  better  (or  at  least  equivalent)  approx- 
imation of  the  actual  harvest  in  western  Montana  than 
did  the  "best"  model  for  two-thirds  of  the  study  years. 
Since  a  separate  model  for  eastern  Montana  was  not 
developed,  we  could  not  assess  whether  the  sum  of 
eastern  plus  western  Montana  provided  a  better  approx- 
imation of  the  Montana  harvest  than  did  the  Montana 
"best"  model.  We  expect,  however,  that  this  would  be 
the  case. 

Finally,  we  present  results  of  analyses  that  modeled 
each  area's  timber  cut  by  the  variables  contained  in  the 
"best"  model  for  each  other  area.  This  was  done  to 
determine  the  extent  to  which  the  need  (perceived  or  re- 
quired) for  model  consistency  between  regions  entailed  a 
loss  of  predictive  ability,  measured  by  R2.  As  planning 
efforts  in  forestry  expand,  there  is  a  tendency  to  require 
procedures  that  ensure  uniformity  and  comparability. 


Results  shown  in  table  4  indicate  a  mixed  response 
again.  The  R2  values  shown  are  frequently  exaggerated 
since  they  correspond  to  models  that  may  contain 
statistically  insignificant  variables.  The  northwestern, 
southwestern,  and  western  Montana  models  all  explained 
an  average  of  72  percent  of  the  variation  in  harvest 
while  the  Montana  model  averaged  about  10  percent 
less.  However,  although  the  northwestern  and  western 
Montana  models  produced  results  that  varied  only  by 
about  22  percent  between  regions,  the  Montana  model 
produced  results  that  varied  by  about  37  percent  and 
the  southwestern  model  by  43  percent.  These  results 
suggest  that  consistency  can  be  achieved  only  at  the 
expense  of  predictive  power,  a  substantial  loss  in  some 
cases.  The  southwestern,  northwestern,  and  western 
Montana  models  all  generated  an  average  loss  of  about 
14  percent  per  (outside)  region.  The  Montana  model  has 
the  largest  loss,  averaging  about  23  percent  per  region. 
Regardless,  either  the  largest  or  smallest  decrease  would 
generally  be  considered  to  be  an  appreciable  loss  in 
predictive  power. 

The  consequence  of  model  consistency  can  be  shown 
by  subsequent  analysis  of  the  southwestern  model  when 
applied  to  other  regions.  Three  sets  of  regression  models 
are  shown,  one  for  each  region  using  the  variables  in  the 
southwestern  model: 

HtMT=      59.6  +  4.0(P)t  -  62.7(Res)t  +  9.4(Unemp)t  .  .  .  R2  =  70.5 
(2.90,  .05)       (n.s.)  (n.s.) 
99.9  +  4.4(P)t  -  56.3  (Res)t  .  .  .  R2  =  68.1 

(4.64,  .01)  (n.s.) 

51.9  +  4.8(P)t  .  .  .  R2  =  67.8 

(5.42,  .01) 


HWM. 


112.1  +  1.9(P)t  -  142.4(Res)t  +  5.7(Unemp)t.  .  .  R2  =  80.2 
(3.72,  .01)  (-1.73,  =.10)  (n.s.) 


HNW= 


78.4  +  2.4(P)t  -  59.3(Res)t 
(5.98,  .01)  (n.s.) 

33.0  +  2.4(P)t 

(6.52,  .01) 

83.0  +  0.8(P)t  -  56.8(Res)t  +  0.7(Unemp)t 
(n.s.)  (n.s.)  (n.s.) 

81.2  +  0.9(P)t  -  51.0(Res)t 
(2.52,  .05)  (n.s.) 

32.3  +  l.l(P)t 

(3.07,  .01) 


.  .  R2  =  76.9 
.  .  R2  =  75.2 
.  .  R2  =  46.6 
.  .  R2  =  46.4 
.  .  R2  =  40.2 


Table  4.— Comparison  of  explanatory  ability  (R2)  of  "best"  model  for  each  region  applied  to  other  regions 


Region 
"best"  model 

"Best" 

model  applied  to 

for 

Northwestern 

Southwestern 

Western 

Montana 

Average 

Montana 

43.5 

62.2 

75.3 

80.7 

65.4 

Western 

76.3 

61.2 

83.7 

67.9 

72.3 

Southwestern 

46.6 

89.3 

80.2* 

70.5 

71.7 

Northwestern 

76.5 

61.3 

83.5 

67.9 

72.3 

5 


All  variables,  t-statistics,  and  significance  levels  are  as 
given  earlier.  These  analyses  rather  clearly  show  that  if 
statistical  significance  (at  a  =  0.10)  of  regression  coeffi- 
cients is  a  requirement  in  model  selection,  the  three- 
variable  southwestern  model  is  not  feasible  in  any  case, 
other  than  the  southwest.  In  conjunction  with  the  other 
variables,  the  unemployment  rate  variable  is  not  signifi- 
cant [(n.s.)]  in  all  cases.  When  the  data  are  fit  to  the  two- 
variable,  (P  and  Res)  model,  the  residual  (Res)  variable  is 
again  not  significant  in  any  of  the  three  regions.  Only 
the  price  variable  (P)  is  statistically  significant  in  all 
regions.  Stated  differently,  starting  from  the  three- 
variable  southwestern  model,  the  only  model  that  could 
be  consistently  applied  to  all  regions  (with  the  signifi- 
cant coefficient  requirement)  is  the  single-variable  price 
model.  In  fact,  that  model  is  the  best,  consistent,  single- 
variable  model  in  light  of  all  possible  variables. 

DISCUSSION 

Three  general  conclusions  derive  from  this  study. 
First,  highly  satisfactory  models  can  be  developed  to 
predict  NIPF  timber  harvest,  if  the  appropriate 
variables  have  been  quantified.  In  particular,  results 
would  have  been  much  less  satisfactory  without  use  of 
the  "cut"  version  of  National  Forest  stumpage  price,  ex- 
pressed in  real  dollars.  Without  the  detailed  county-level 
records  on  NIPF  timber  harvest,  our  analysis  of  Mon- 
tana regions  simply  could  not  have  taken  place.  This 
type  of  record  is  not  consistently  available  in  all  States. 
Second,  the  ability  to  develop  better  predictive  models 
for  small  areas  compared  to  larger  areas  was  not 
uniform.  Sometimes  it  worked;  sometimes  it  did  not. 
This  suggests  that  the  assumed  advantage  of  regional 
forecasting,  reducing  variability,  does  not  always  occur. 
Third,  use  of  the  same  model  to  predict  timber  harvest 
for  different  regions  has  both  advantages  and  disadvan- 


tages. The  clear  advantages  are  consistency  and  the  at- 
tendant ease  in  understanding  the  model  and  inter- 
preting results.  But  these  advantages  can  only  be 
secured  by  a  loss  of  predictive  ability.  Results  presented 
indicate  that  the  loss  can  be  substantial  or  nearly 
inconsequental. 

Throughout  this  report,  the  coefficient  of  multiple 
determination  (R2)  has  been  used  as  the  basis  for  com- 
parisons and  for  judging  model  superiority.  While  a  con- 
venient measure,  R2  alone  is  probably  an  inadequate 
basis  for  such  judgments.  Other  measures  such  as  ad- 
justed R2,  standard  errors,  and  ease  of  quantifying  in- 
dependent variables  should  also  be  considered.  Moreover, 
although  some  of  the  R2  values  presented  were  substan- 
tially different  from  others,  some  were  quite  close.  An  R2 
of  80.2  may  be  statistically  different  from  83.7,  but  the 
difference  could  be  judged  trivially  small  in  a  practical 
application.  The  ultimate  determination  that  one  model 
is  "close  enough"  to  another  is  a  judgment  that  must  be 
based  on  the  needs  and  priorities  of  the  user. 
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